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Abstract . — A parsimony analysis on morphological characters was per- 
formed to estimate the phylogenetic relationships of the taxa of Pinus subsect. 
Australes. The Adams consensus tree placed the Caribbean species as a mono- 
phyletic clade with P. rigida—P. serotina as its sister taxon. Based on this 
phylogeny, area cladograms were constructed and compared to geologic clado- 
grams constructed from plate-tectonic evidence. This comparison and an an- 
cestral area analysis indicate that colonization of the tropics most likely oc- 
curred from Florida to Hispaniola, rather than by the circumferential-Gulf 
route. Subsequent dispersal events to Central America, Cuba, and the Bahamas 
are proposed to explain the geographic distribution of P. caribaea. Ecological 
comparisons within subsect. Australes found that sister species are not syntopic 
and that syntopic species are not sister species. Although some North American 
sister species are ecologically quite different, there is low ecological diversity 
among the Caribbean species. 



Pinus is one of the most widespread gen- 
era of plants in the northern hemisphere 
(Mirov 1967:307-308; Strauss & Doerksen 
1990). Its species are found from Central 
America and Sumatra to the Arctic circle 
(Little & Critchfield 1969, Strauss & 
Doerksen 1990, Farjon 1996), occupy xeric 
to mesic habitats, and comprise one of the 
dominant vegetation types on the earth. 
They are important economically, being 
used for fuel and lumber, and are frequently 
planted for commercial purposes in parts of 
the world where they do not naturally occur 
(Mirov 1967:451). While often considered 
early successional species well-adapted to 
poor soils (Govindaraju 1984), pines are ac- 
tually ecologically diverse and are found 
from sea level to the timberline, from sea- 
sonally wet savannas to deserts, and from 
monotypic stands to multispecies climax 
forests where there is co-occurrence with 



hardwoods. McCune (1988) characterized 
the ecological diversity for North American 
pines and defined five ecological groups: 
(1) fire-resistant species, (2) mesophytic 
shade-tolerant species, (3) stress-tolerant 
species, (4) fire-resilient species, and (5) 
southern mesic species. 

Other studies provide a more quantita- 
tive assessment of the interaction between 
pines and their environment. In particular, 
studies of bark ontogeny and fire have 
shown most pines to be particularly well 
adapted to surviving surface fires (Har- 
mon 1984, D. C. Adams 1994). Recent 
work by Jackson and Adams demonstrates 
the evolution of negative bark allometry 
(D. C. Adams & Jackson 1995) in species 
whose habitats are characterized by fre- 
quent surface fires. This finding is consis- 
tent with models of defensive structure 
evolution, where resources are allocated 
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earlier to defense as the likelihood of a 
mortality factor increases. 

Although the ecological roles of pines 
are relatively well defined, aspects of the 
classification and systematics of the genus 
Pinus have been uncertain (Mirov 1967: 
540, Strauss & Doerksen 1990). Shaw (1914) 
published the first major taxonomic work 
on pines, where he recognized two subgen- 
era, Strobus and Pinus. Pilger (1926) pre- 
sented a second classification of the pines, 
classifying many of Shaw’s (1914) varieties 
as species and further dividing the subgen- 
era into eleven sections. Because of his 
heavy reliance on needle number, however, 
Pilger’s classification is considered to be a 
step backwards (Mirov 1967:526). Duffield 
(1952) revised the classification of the 
pines, using information from hybridization 
studies. He divided Shaw’s “Group Austra- 
les ” into two “Groups,’’ the eastern species 
(XI) and the western species (XII). In 
Shaw’s (1914) classification, most of the 
southern pines (P. echinata, P. elliottii, P. 
glabra, P. palustris, and P. taeda) as well 
as the Caribbean and many western species, 
were included in the “Group Australes" . 
Duffield also added three species from 
Shaw’s “Group Insignes ” to his “Group 
XI”: P. rigida, P. serotina, and P. pungens. 
Little & Critchfield (1969) reviewed the 
many classification schemes for Pinus. In 
particular, they restored the name Australes, 
at the subsectional rank, to Duffield’s 
“Group XI,” which now included the eight 
southern pines and three Caribbean species. 

With the advent of phylogenetic meth- 
odology, many subsections of the genus Pi- 
nus have been reexamined in an evolution- 
ary context, so that the relationships among 
many of the North American taxa are now 
better known (Wheeler et al. 1983, Strauss 
& Doerksen 1990, Malusa 1992, Govinda- 
raju et al. 1992, Krupkin et al. 1996). The 
species comprising subsect. Australes, how- 
ever, have not been examined in a phylo- 
genetic context. In particular, it is of interest 
to determine whether the Caribbean taxa 
comprise a monophyletic clade and, if so. 



what its relationship is to the North Amer- 
ican species. Mirov (1967:555) thought that 
insular varieties of P. caribaea originated 
from a Central American progenitor, but 
without a phylogenetic analysis he could 
not specify whether this was the original 
colonization of the Caribbean islands by 
subsect. Australes or was part of a subse- 
quent radiation within the Caribbean re- 
gion. Farjon (1996) on the other hand, sug- 
gested that the Caribbean taxa were of a 
North American origin, and that the Central 
American P. caribaea originated from Ca- 
ribbean immigrants. Because of these out- 
standing questions, we felt that an exami- 
nation of the relationships of the taxa in 
subsect. Australes was needed, and there- 
fore performed a phylogenetic analysis of 
the southern pines using morphological 
characters. In addition, we generated an 
area cladogram for the geographic regions 
occupied by the taxa and compared this to 
geologic evidence in order to better under- 
stand the vicariance and dispersal of the 
taxa in the Caribbean. 

Methods 

Phylogenetic analysis. — The eight spe- 
cies of subsect. Australes in the eastern 
United States and certain geographically 
defined OTUs of the three Caribbean spe- 
cies were used in this study (Appendix I). 
Pinus caribaea from the Bahamas, Cuba, 
Belize, and Honduras-Nicaragua were treat- 
ed as separate OTUs. Because P. caribaea 
in the uplands of Belize and Honduras may 
be subject to hybridization with P. oocarpa 
(Williams 1955) we used only specimens 
from lowland localities. Pinus elliottii var. 
densa was also treated as a separate OTU. 
The subsection Australes is included within 
the hard pines (subg. Pinus ) as classified by 
Little & Critchfield (1969). While the re- 
lationships among the hard pines are still 
somewhat controversial (see Strauss & 
Doerksen 1990, Govindaraju et al. 1992), 
we chose to use two taxa (P. virginiana and 
P. clausa) from subsect. Contortae as out- 
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Table 1. — Data matrix for the taxa of Pinus used in 
this study. Descriptions of characters and character 
codes are found in Appendix II. 



Species 


Character states 


P. caribaea — Bahamas 


0000101011 


1101311112 7 


P. caribaea — Belize 


0001111011 


1101411111 8 


P. caribaea — Cuba 


0000111001 


1100511111 9 


P. caribaea — Hondu- 






ras — Nicaragua 


0001111011 


1101411112 6 


P. clausa 


0010101110 3100100000 2 


P. cubensis 


00001 1 101 1 


1111211014 7 


P. echinata 


0000011111 


1100200000 2 


P. elliottii var. elliottii 


1000110011 


1100211012 2 


P. elliottii var. densa 


1000111011 


1100211013 3 


P. glabra 


0000011111 


1100100000 2 


P. occidentalis 


0000101011 


1101611113 3 


P. palustris 


1001100000 2000401113 0 


P. pungens 


1101100000 3100201005 0 


P. rigida 


0000111011 


3100411102 1 


P. serotina 


0000111011 2101511102 0 


P. taeda 


1101100000 3100301100 4 


P. virginiana 


0010101110 3100101000 5 



groups, based on the phylogenies of Gov- 
indaraju et al. (1992) and Farjon (1996). 

Twenty one morphological characters 
were scored from herbarium specimens 
(LAF, MO, NY), with supplementary infor- 
mation from North American specimens 
collected by us and from literature sources 
(Faijon & Styles 1997, Radford et al. 
1964); species were assigned values based 
on an average of several specimens. Four- 
teen characters were based on cone mor- 
phology and seven on needle morphology 
(Appendix II). Four of these characters 
were coded as ordered, multi-state charac- 
ters (Table 1; Appendix II). Binary scoring 
was employed for the remainder. While rec- 
ognizing that some of the characters are 
quantitative (Stevens 1996), we believe bi- 
nary scoring adequately approximates char- 
acter states among taxa where we have used 
it. We generated a phylogenetic hypothesis 
through Wagner parsimony using the 
branch and bound algorithm in PAUP ver- 
sion 3.1 (S wofford 1991). Wagner parsi- 
mony attempts to reconstruct an evolution- 
ary tree by minimizing the number of 
changes of the character states along the 



tree (Kluge & Farris 1969, Strauss & 
Doerksen 1990). If one assumes that a spe- 
cies from a closely related group reflects 
previous character states, it can be used as 
an outgroup to root the tree and to polarize 
the character states as well (Wiley 1981). 

Biogeographic analysis . — Comparing 
the current distributions of taxa with their 
phylogenetic history to elucidate patterns, 
as well as to evaluate the relative plausibil- 
ities of vicariance and dispersal, is the pur- 
pose of cladistic biogeography (Nelson & 
Platnick 1981, Wiley 1988, Morrone & Cri- 
si 1995). To examine the biogeography of 
subsect. Australes in a phylogenetic con- 
text, we used the method outlined by 
Brooks & McLennan (1991). After the phy- 
logenetic relationships were estimated, the 
geographic areas containing the taxa were 
defined: eastern North America (NA), the 
Bahamas (B), Cuba (Cb), Hispaniola (H), 
the Yucatan peninsula (Y), and Honduras- 
Nicaragua (Ch). The relationships of these 
geographic areas were then determined by 
representing the relationships of the taxa by 
a matrix containing additive binary codes, 
replacing the species with their respective 
geographic areas, and performing a parsi- 
mony analysis. Brooks Parsimony Analysis 
(BPA) has been criticized for its treatment 
of widespread taxa as synapomorphous for 
areas (Kluge 1988) and as a method for 
generating general area cladograms (Nelson 
& Ladiges 1991). Our goal was to evaluate 
alternative historical biogeographies rather 
than to produce a general area cladogram. 
Having treated the species population of 
each major Caribbean locale as an OTU, 
our data set did not, by definition, have 
widespread taxa, but there was a redundant 
distribution with two OTUs in Cuba. We 
chose to derive area cladograms by BPA 
with (inclusive ORing: Cressey et al. 1983) 
and without redundant distributions in order 
to evaluate the redundancy (Brooks 1990). 
The ancestral area of the Caribbean taxa 
was estimated by the technique of Bremer 
(1992). This involves comparing areas in 
regard to the numbers of necessary gains 
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and losses in the area cladogram by Camin- 
Sokal parsimony and under the assumption 
that each area is ancestral. The putative an- 
cestral area best supported by the biogeo- 
graphic evidence is that having the highest 
ratio of gains to losses. 

Phylogenetically-based cladograms are 
intended for comparison with geological 
area cladograms. Much work has been done 
on the biogeography of the Caribbean flora 
and fauna (e.g., Rosen 1976, 1978, 1985; 
Guyer & Savage 1986; Page & Lydeard 
1994; Hedges et al. 1994). However, due to 
the complex geologic history of the region 
(see Pindell & Barrett 1990), many of the 
interpretations of biogeography are highly 
controversial and out of date. In light of the 
increased knowledge of Caribbean geology 
(Burke 1988, Pindell & Barrett 1990), we 
felt that using geologic cladograms from 
previous biogeographic work might be un- 
wise. We therefore compared the area 
cladograms for subsect. Australes to geo- 
logic cladograms generated from recent tec- 
tonic evidence on the relationships of the 
geologic regions of the Caribbean. 

Results 

Phylogenetic analysis. — We found three 
most parsimonious trees from our cladistic 
analysis, each containing 81 steps and a 
consistency index (Cl) of 0.469 (Fig. 1). 
When uninformative characters were ex- 
cluded, the consistency index was reduced 
by only 0.010 to Cl = 0.459, implying that 
most characters were phylogenetically in- 
formative. From these three most parsimo- 
nious trees we generated a consensus tree, 
using the procedure described by E. N. Ad- 
ams (1972, 1986). This tree contains only 
that information present in all rival tree to- 
pologies, and is thus a conservative esti- 
mate of the true topology. 

Based on this phylogeny, subsect. Aus- 
trales is divided into several distinct sub- 
clades. The smaller subclade is a polytomy 
containing P. taeda, P. pungens, and P. pa- 
lustris. A second polytomy is located deep- 



er in the phylogeny and contains the smaller 
subclade, the larger subclade, and P. echin- 
ata and P. glabra as single taxa. The re- 
maining taxa are all found within the larger 
subclade. That the phylogeny locates P. rig- 
ida and P. serotina as sister taxa accords 
with the proposal by Smouse and Saylor 
(1973), who consider them conspecific. The 
most significant aspect of this phylogeny is 
that the Caribbean species are monophyletic 
within the larger subclade. 

Within the larger subclade, it is notewor- 
thy that P. elliottii var. densa is the sister 
taxon to P. elliottii var. elliottii. Based on 
the chemical composition of its turpentine, 
Mirov et al. (1965) and Mirov (1967:555) 
proposed that P. elliottii var. densa was 
more closely related to the Caribbean spe- 
cies, and was a recent arrival to Florida. 
Our findings suggest that it is in fact more 
closely related to the mainland slash pine 
than to any Caribbean taxon. In addition, 
our phylogeny suggests that P. cubensis is 
more closely related to P. caribaea than to 
P. occidentalis, in contradiction to the pro- 
posal by Mirov (1967:232). Our findings 
place P. occidentalis basal to the Caribbean 
subclade. 

Biogeographic analysis. — Much of the 
history of the Caribbean region is still un- 
known, as is evident from the numerous 
tectonic models of the region (Perfit & Wil- 
liams 1989). However, the present study fo- 
cuses on taxa in the Greater Antilles, the 
Bahamas, the Yucatan peninsula, and Hon- 
duras-Nicaragua, where there seems to be a 
general consensus of opinion concerning 
the geologic history. We therefore present a 
brief account of the geologic history of the 
region (summarized in Fig. 2) and use this 
to generate geologic cladograms for the 
taxa. 

Mexico collided with the North Ameri- 
can plate during the Jurassic (Burke 1988, 
Pindell & Barrett 1990). The Greater An- 
tilles were part of a larger body called the 
Great Arc, which originated in the Pacific 
during the late Cretaceous and migrated 
northeast into the Atlantic in the Paleocene 



VOLUME 1 10, NUMBER 4 



685 




P. taeda 
P. pungens 
P. palustris 
P. elliottii 
P. elliottii-densa 
P. rigida 
P. serotina 

P. occidentalis-Hispaniol 
P. cubensis 
P. caribaea-Bahamas 
P. caribaea-Cuba 
P. caribaea-Belize 
P. caribaea-Honduras 
P. echinata 
P. glabra 
P. virginiana 
P. clausa 



Fig. 1 . Adams consensus tree found from three most parsimonious trees based on twenty-one morphological 
characters of Pinus spp. 



(Burke 1988, Pindell & Barrett 1990). The 
Great Arc broke into three segments some- 
time in the late Cretaceous to early Paleo- 
cene (Burke 1988). The northern segment 
became the Greater Antilles, the central be- 
came the Lesser Antilles, and the southern 
segment is presumed to have collided with 
the South American continent. 

During the late Paleocene, the Greater 
Antilles began to collide with the Bahaman 
plate (Pindell & Dewey 1982, Burke 1988, 
Pindell & Barrett 1990). However, Burke 
(1988) claims that prior to this, the Greater 
Antilles collided with the Yucatan. Others 



(Pindell & Barrett 1990) claim that al- 
though there is geologic evidence for such 
a collision, it cannot be determined whether 
it was the Greater Antilles or some other 
geologic body that collided with the Yuca- 
tan. Strike-slip faults began to separate 
Cuba and Hispaniola in the mid-Eocene 
(Pindell & Barrett 1990), and the Honduras- 
Nicaraguan block (called Chords) began to 
collide with the Yucatan shortly after in the 
Miocene (Perfit & Williams 1989). Though 
Pindell and Barrett place this collision 
slightly earlier in the Oligocene, this does 
not affect our geologic cladograms. Based 
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Fig. 2. Tertiary paleogeography of the Caribbean (redrawn from Perfit and Williams 1989). A. Late Creta- 
ceous. B. Paleocene. C. Mid-Eocene. D. Late Miocene-Pliocene. AR = Aves Ridge, B = Bahaman Platform, 
CA = Cuban Arc, CH = Chords, H = Hispaniola, J = Jamaica, L = Lesser Antilles Arc, NH = North Island 
of Hispaniola, SH = South Island of Hispaniola. 



on this geologic information, we have con- 
structed two area cladograms (Figs. 3a, b). 
The first (Fig. 3a) represents the relation- 
ships of the areas if the Yucatan had not 
collided with the Greater Antilles ( sensu 
Pindell & Barrett 1990), and the second 
(Fig. 3b) if it had {sensu Burke 1988). 

The taxon-based area cladogram derived 
with inclusive ORing has Hispaniola as the 
sister area to the other Caribbean areas and 
situates Cuba as the sister area to the Yu- 
catan, in a derived position (Fig. 4a). The 
taxon-based area cladogram that treated 



Cuba as two separate areas historically (ex- 
clusive ORing) has a pattern of branching 
identical to that of the taxon cladogram, but 
with Cuba as a sister area to both the Yu- 
catan and the Bahamas (Fig. 4b). The an- 
cestral area analysis (AAA) of Bremer 
(1992) was carried out in two ways: on an 
area cladogram with terminal sister taxa not 
grouped, which corresponds to Fig. 4b; and 
on an area cladogram that considered there 
to be no node separating the two OTUs of 
P. caribaea of Cuba and Belize. The logic 
of combining the areas of the terminal sister 
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A 



NA Y Ch B Cb H 





Fig. 3. Geologic cladograms based on tectonic 
evolution in the Caribbean. (A) corresponds to hy- 
pothesis of Pindell and Barrett (1990), where the Yu- 
catan did not collide with the Greater Antilles; (B) 
corresponds to Burke’s (1988) hypothesis, which in- 
cludes this collision. 

taxa from this node is that the taxa are con- 
specific. For both approaches to AAA, His- 
paniola is the most probable ancestral area, 
and Cuba is the next most probable (Table 
2 ). 

The positions of Yucatan-Chortis and 
Hispaniola differ between the geological 
and taxon-based area cladograms (Figs. 3 
and 4). In the geological area cladograms, 
Yucatan-Chortis is either the sister area to 
North America or the basal area in the Ca- 
ribbean, whereas the taxon-based area 
cladograms have Yucatan-Chortis separated 
and both in more derived positions. Hispan- 
iola is placed near the origin of the Carib- 
bean subclade by the taxon-based area 
cladogram, but in a derived position by the 
geological data. We conclude that the 
cladogram based on geological contact and 
separation corresponds poorly to the se- 
quence of areas colonized by the Caribbean 
taxa of subsect. Australes. 

The derived placement of Cuba in the 



A 



NA Ch B Y Cb H 




NA Ch B Cb2 Y Cbl H 




Fig. 4. Area cladograms of six geographic locali- 
ties based on the phylogeny for Pinus subsect. Austra- 
les using: (A) inclusive ORing, and (B) exclusive OR- 
ing. Letters correspond to the different geographic 
regions: NA = southern United States, Ch = Chords, 
B = Bahamas, Y — Yucatan, Cb — Cuba, H = His- 
paniola. 

taxon-based area cladogram calculated by 
inclusive ORing is due to there being two 
derived taxa ( P . cubensis and P. caribaea ) 
on Cuba that determine the location of 
Cuba on the area cladogram (Brooks & 

Table 2. — Estimation of ancestral area for the Ca- 
ribbean subclade of the subsection Australes . Values 
not in parentheses are for cladogram that reduces con- 
specific sister taxa. Values in parentheses are for clado- 
gram without reduction of conspecific sister taxa. G = 
number of necessary gains under forward Camin-Sokal 
parsimony. L = number of necessary losses under re- 
verse Camin-Sokal parsimony. AA = G/L rescaled by 
division by the largest G/L. AA = mean of AA’s cal- 
culated by the two reduction alternatives. 



Area 


G L G/L AA AA 


Hispaniola 

Cuba 

Chords 

Yucatan 

Bahamas 


1 (1) 1 (1) 1.00 (1.00) 1.00 (1.00) 1.00 

2 (2) 3 (4) 0.67 (0.50) 0.67 (0.50) 0.59 

1 (1) 2 (2) 0.50 (0.50) 0.50 (0.50) 0.50 

1 (1) 3 (4) 0.33 (0.25) 0.33 (0.25) 0.29 

1 (1) 4 (4) 0.25 (0.25) 0.25 (0.25) 0.25 
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McLennan 1991:212). The presence of de- 
rived taxa and the lack of ancestral taxa on 
Cuba can be explained through intra-Cuban 
differentiation combined with colonization 
of other Caribbean areas from Cuba. The 
AAA clearly supports this explanation. 

Discussion 

Based on our phylogenetic hypothesis of 
the southern pines, several ecoevolutionary 
trends can be recognized. First, sister spe- 
cies are not syntopic. For example, the sis- 
ter clades of P. rigida and P. serotina as 
well as P. cubensis and P. caribaea repre- 
sent taxa that are allopatric. Pinus pungens 
and P. palustris are allopatric as well and 
occupy different habitats. Second, species 
that are syntopic are not sister species. Ex- 
amples of this are P. palustris and P. el- 
liottii in the eastern Gulf coastal plain flat- 
woods, P. palustris and P. echinata in the 
uplands of the western Gulf coastal plain, 
and P. taeda and P. echinata in the south- 
ern Piedmont. 

Major ecological divergence is found be- 
tween some sister species. For example, P. 
pungens grows in dense monospecific 
stands in xeric montane locations and main- 
tains dominance by mass recruitment after 
stand-replacing crown fires; its bark grows 
with positive allometry (Adams 1994). In 
comparison, P. palustris grows as widely 
spaced individuals in parkland vegetation 
subject to frequent surface fires, and early 
life history stages of this species are fire- 
resistant (McCune 1988), as typified by its 
“grass” stage and negative allometry of 
bark. Curiously, there is a general lack of 
ecological divergence in the Caribbean sub- 
clade of subsect. Australes (Smith 1954). 
This subclade is less diverse in that there 
are no mountain ridge specialists like P. 
pungens, there are no shade-tolerant species 
like P. glabra, and there are no species with 
a fire-resistant seedling stage like P. palus- 
tris or P. elliottii var. densa. 

The purpose of cladistic biogeography is 
to attempt to discern whether current dis- 



tribution patterns are the result of vicariance 
or dispersal. In subsect. Australes this issue 
is related to whether the initial colonization 
of Caribbean islands was from eastern 
North America or from Central America. 
Both the ancestral area analysis and the fact 
that neither Yucatan nor Chortis branch low 
on the area cladogram suggest that the 
Greater Antilles were the first sites of col- 
onization in the Caribbean. This could have 
occurred by dispersal from North America 
in the late-Tertiary or Quaternary, or by vi- 
cariance from Central America if a com- 
ponent of Hispaniola collided with, and 
then separated from, Yucatan during the Pa- 
leocene (Burke 1988). The issue is one of 
timing and therefore cannot be decided con- 
clusively without molecular data on the di- 
vergence time between North American and 
Caribbean taxa. However, several facts ar- 
gue in favor of dispersal to Hispaniola. For 
example, there is no evidence in Central 
America of a Caribbean ancestor; P. cari- 
baea cannot play this role given its derived 
status. The Caribbean subclade was derived 
from an ancestor that shared many traits, 
particularly of cones, with P. serotina. Cur- 
rently P. serotina has an Atlantic-East Gulf 
Coastal Plain distribution, not being found 
west of Mobile Bay. This puts P. serotina 
geographically closer to the late-Tertiary 
Greater Antilles than to the putative contact 
between the Yucatan and the proto-Greater 
Antilles. In addition, the lack of ecological 
diversification in the Caribbean clade could 
suggest an occupation of the region too re- 
cent to have occurred via Paleocene vicar- 
iance. Farjon (1996) found that P. occiden- 
talis and P. caribaea var. hondurensis 
formed a sister group to all the other neo- 
tropical pines he considered, a result con- 
sistent with derivation from a southeastern 
North American ancestor. 

Dispersal of Australes to the Greater An- 
tilles may have been part of a pattern of 
colonization by xeric-adapted biota. R. P. 
Adams (1989) presented evidence from leaf 
morphology and volatile leaf oils that 
strongly- indicates an origin of West Indian 
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Juniperus species from eastern North 
America, rather than from Mexico. Buck 
(1990) suggested that several xeric bryo- 
phyte species found in North America and 
upland Hispaniola colonized the island via 
dispersal during glacial episodes of the 
Pleistocene when over-water distances were 
less due to lower sea levels (Gascoyne et 
al. 1979) and when savanna habitats were 
more widespread than at present because of 
a more xeric climate (Pregill & Olson 
1981). Such a scenario could apply to col- 
onization of the Greater Antilles by subsect. 
Australes as well as to inter-island move- 
ment and colonization of Central America. 

Mirov (1967:555) hypothesized that P. 
caribaea colonized the Caribbean islands 
from Central America. The area cladogram, 
with Chortis as the sister area to all other 
areas occupied by P. caribaea, supports this 
direction of colonization in the species. Dis- 
persal is implicated because no connection 
is known to have existed between Chortis 
and the Greater Antilles. The only possibil- 
ity of vicariance within the P. caribaea 
clade would be between the taxonomic sub- 
units in Cuba and Yucatan. However, any 
such vicariance would date from early Ter- 
tiary (Burke 1988), and this appears too 
longstanding for the small level of differ- 
entiation between the taxa. It is difficult to 
ascribe any of the divergences within the 
Caribbean subclade to vicariance. 

It was believed by Mirov et al. (1965) 
that P. elliottii var. densa was a recent im- 
migrant to Florida from the Caribbean and 
that it was therefore more closely related to 
the Caribbean species. Our results place P. 
elliottii var. densa as the sister taxon to P. 
elliottii var. elliottii, suggesting that it is 
more closely related to the mainland spe- 
cies. Squillace (1966) demonstrated multi- 
character clinal variation between P. elliot- 
tii var. elliottii and P. elliottii var. densa in 
central Florida. We also found that P. cub- 
ensis is more closely related to P. caribaea 
than to P. occidentalis. Mirov (1967:232) 
had proposed that P. cubensis and P. occi- 
dentalis were closely related taxa. Clearly 



more work is needed to determine their re- 
lationship. 

Klaus (1980) suggested that, for compar- 
isons through nodes deep within the phy- 
logeny of Pinus, cones display primitive 
character states at the apex and derived 
states at the base. Consideration of intra- 
cone character variation in the context of 
the subsect. Australes cladogram provides 
evidence of this phenomenon at a finer phy- 
logenetic scale. Pinus glabra, P. elliottii, 
and P. rigida often have flattened, distally- 
pointed umbos at the cone base and more 
erect umbos at the apex. Within the clado- 
gram, erect umbos characterize the out- 
group and the P. taeda-P. pungens-P. pa- 
lustris subclade, whereas flattened umbos 
are typical in P. serotina and the Caribbean 
subclade. In P. caribaea var. hondurensis, 
the umbo spine tends to be unconnected to 
the keel on basal scales, a character state 
more widespread in P. cubensis cones, but 
connected to the keel on apical scales, a 
state found throughout the cone in taxa de- 
rived from nodes below P. caribaea var. 
hondurensis. A counterexample, however, 
exists in the relatively erect umbos apically 
in some P. caribaea var. hondurensis. This 
character state is derived if the Caribbean 
clade arose through an ancestor like P. oc- 
cidentalis, yet it appears at the apex. 

Our phylogenetic hypothesis of subsect. 
Australes is a first step in understanding the 
history of the pines of the southern United 
States and the Caribbean. Based on this 
phylogeny, we have identified possible dis- 
persal and vicariance events and have pro- 
vided a framework on which future studies 
may be based. While it is appealing to in- 
terpret phylogenetic hypotheses as recon- 
structions of evolutionary history, it must 
be stressed that they are only estimations of 
the true topology based on the available 
data. Studies have shown that most phylo- 
genetic methods can perform rather poorly 
in their estimation of true tree topology (Fi- 
ala & Sokal 1985; Rohlf et al. 1990). We 
therefore present this phylogenetic hypoth- 
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esis as a tentative estimation of the rela- 
tionships of the southern pines. 
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Appendix I 

Herbarium specimens from which 
characters were scored. 

Pinus caribaea Morelet var. bahamensis Barrett and 
Golfari. BAHAMAS. Grand Caicos: Correll 49463 
(NY). New Providence: Degener 18753 (NY). P. car- 
ibaea var. caribaea Barrett & Golfari. CUBA. Pinar 
del Rio: Jack 8673 (NY), Leon & Charles 4935 (NY). 
P. caribaea var. hondurensis B. & G. BELIZE. Belize: 
Croat 24005 (MO), Kluge s.n. (MO). Corozal: Crane 
315 (MO). Orange Walk: Lundell 677 (MO). Stann 
Creek: Stevenson 1 128 (MO). Toledo: Gentle 3690 
(MO). HONDURAS. Gracias a Dios: Clewell 4511 
(MO), Nelson & Hernandez 1004 (MO). NICARA- 
GUA. Zelaya: Marshall & Neill 6559 (MO), Seymour 
3650 (MO), Stevens 7638 (MO), Stevens 7753 (MO), 
Stevens 21679 (MO), Vincelli 625 (MO). P. clausa 
(Chapm.) Vasey. USA. Florida: Franklin Co., Godfrey 
69282 (LAF); Polk Co., Shuey 2317 (LAF); Walton 
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Co., Smith 2071 (LAF). P. cubensis Griseb. CUBA. 
Oriente: Ekman 3829 (NY), Shafer 4174 (NY). P. 
echinata Mill. USA. Louisiana: Bienville Par., Thieret 
16836 (LAF), Westling 141 (LAF); Evangeline Par., 
Reese & Reese 1537 (LAF), Thieret 22256 (LAF); St. 
Helena Par., Allen 1499 (LAF). P. elliottii Engelm. var. 
densa Little & Dorman. USA. Florida: Dade Co., 
Small & Carter 1249 (NY); Monroe Co., Brizicky & 
Stem 378 (NY). P . elliottii var. elliottii L. & D. USA. 
Louisiana: Lafayette Par., Thieret 17423 (LAF); St. 
Tammany Par., Thieret 16770 (LAF), Thieret 21879 
(LAF); Tangipahoa Par., Thieret 16758 (LAF). P . gla- 
bra Walt. USA. Louisiana: Livingstone Par., Thieret 
16743 (LAF); St. Helena Par., Allen 1500 (LAF), St. 
Tammany Par., Lynch 940 (LAF); Washington Par., 
Thieret 16782 (LAF). P. occidentalis Swartz. DOMIN- 
ICAN REPUBLIC. Zanoni, Mejia, Pimentel & Garcia 
32443 (NY). HAITI. Nash 825 (NY). P. palustris Mill. 
USA. Louisiana: Allen Par., Thieret 10389 (LAF);- 
Beauregard Par., Thieret 16968 (LAF); Nachitoches 
Par., Thieret 16811 (LAF), Thieret 17035 (LAF); St. 
Helena Par., Thieret 17298 (LAF); Tangipahoa Par., 
Thieret 16782 (LAF); Winn Par., Thieret 16817 (LAF). 
P. pungens Lamb. USA. Pennsylvania: Franklin Co., 
Adams s.n. (LAF), Adams 25 (LAF). P. rigida Mill. 
CANADA. Quebec: Lemieux 1281 (LAF). USA. 
Maryland: Frederick Co., Windier & Stastny 3900 
(LAF). North Carolina: Buncombe Co., Dunton s.n. 
(LAF); Jackson Co., Duncan 22793 (LAF). Pennsyl- 
vania: Franklin Co., Adams s.n. (LAF). West Virginia: 
Hardy Co., Wratchford s.n. (LAF). P. serotina Michx. 
USA. Alabama: Geneva Co., Krai 33947 (LAF). Geor- 
gia: McIntosh Co., Duncan 20688 (LAF). South Car- 
olina: Charleston Co., Ahles 53156 (LAF), P. taeda L. 
USA. Louisiana: Livingston Par., Thieret 16742 (LAF) 
Thieret 16753 (LAF), Thieret 16754 (LAF); St. Helena 
Par., Thieret 17300 (LAF); Tangipahoa Par., Thieret 
16756 (LAF). P. virginiana Mill. USA. Alabama: 
DeKalb Co., Vincent 1281 (LAF). Georgia: Rabun 
Co., Reade s.n. (LAF). Maryland: Allegany Co., Dun- 
can 22944 (LAF). North Carolina: Orange Co., Ahles 
53117 (LAF). Tennessee: Cheatham Co., Demaree 
49183 (LAF). 

Appendix II 

Twenty-one morphological characters used in this 
study. The first fourteen characters describe mature 



megasporangiate cone morphology and the last seven 

describe needle morphology. 

(1) Sum of maximum cone length and maximum 
cone width: less than 20 cm = 0; greater than 20 
cm = 1. 

(2) Cones subsessile or stalked = 0; sessile = 1. 

(3) Transverse keel or dorsal surface of cone scale 
between apophysis and umbo: depressed at junc- 
tion = 0; continuous at junction = 1. 

(4) Lateral portion of most umbos not strongly ele- 
vated above apophysis = 0; strongly elevated 
above apophysis = 1. 

(5) Umbo: weakly keeled or unkeeled = 0; strongly 
keeled = 1. 

(6) Umbo keel: straight or monotypically curved = 
0; undulating = 1. 

(7) Proximal portion of umbo slopes at: high angle 
or is perpendicular to apophysis surface = 0; low 
angle from spine to proximal edge of umbo = 1. 

(8) Proximal portion of umbo: not concave = 0; con- 
cave = 1. 

(9) Portion of umbo proximal to keel: not larger than 
distal portion = 0; much larger than distal portion 
= 1 . 

(10) Junction of distal margin of umbo and apophysis: 
grooved = 0; not grooved = 1. 

(11) Spine on umbo: unkeeled on apical half = 1; 
keeled on part of apical half but unkeeled at apex 
= 2; keeled on all of apical half = 3. 

(12) Spine on umbo: reflexed = 0; straight or curved 
outward = 1. 

(13) Spine: connected to keel = 0; substantially distal 
to keel = 1. 

(14) Most umbos erect such that apex of spine is not 
pointed distally near the level of the apophyseal 
keel = 0; most umbos flattened and turned dis- 
tally such that apex of spine is pointed distally 
near the level of the apophyseal keel = 1. 

(15) Needle number: 2=1; 2(— 3) = 2; 3(-2) = 3; 
3 = 4; 3(— 4) = 5; 3-5 = 6. 

(16) Hypodermal cells in leaf angles: absent = 0; 
present = 1. 

(17) Hypodermis structure: uniform = 0; biform = 1. 

(18) Endodermal cells: walls normal = 0; walls thick- 
ened = 1. 

(19) Resin canals medial = 0; internal = 1. 

(20) Number of resin canals: coded character. 

(21) Number of stomata: coded character. 



